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Abstract 

We study the evaporation of (4 + n)-dimensional rotating black holes into scalar degrees 
of freedom on the brane. We calculate the corresponding absorption probabilities and 
cross-sections obtaining analytic solutions in the low-energy regime, and compare the 
derived analytic expressions to numerical results, with very good agreement. We then 
consider the high-energy regime, construct an analytic high-energy solution to the scalar- 
field equation by employing a new method, and calculate the absorption probability 
and cross-section for this energy regime, finding again a very good agreement with the 
exact numerical results. We also determine the high-energy asymptotic value of the total 
cross-section, and compare it to the analytic results derived from the application of the 
geometrical optics limit. 



1 Introduction 



Among the motivations for consideration of higher-dimensional theories [U |2] is that the 
leading candidates (String Theory and variants of it), for the unification of gravity with 
the rest of the fundamental interactions at the quantum level, are all formulated in a 
higher- dimensional context. In these models, gravity propagates in D = 4 + n dimensions 
(Bulk), while matter degrees of freedom are confined to live on a 4-dimensional hyper- 
surface (Brane). In models with large extra dimensions [T], the traditional Planck scale 
Mpi ~ 10 18 GeV is only an effective scale, related to the fundamental higher-dimensional 
gravity scale M* through the relation M Pl ~ M™ +2 i? n , where R ~ (Ki) 1 /™ is the effective 
size of the n extra spatial dimensions. If R 3> £pi ~ 10~ 33 cm, the scale M* can be sub- 
stantially lower than Mpi. In that case, trans-planckian particle collisions could probe the 
strong-gravity regime and possibly produce higher-dimensional black holes [3] centered 
at the brane and extending in the bulk. For all the classical laws of black-hole physics to 
still hold, the mass of the black hole Mbh would have to be larger than M*. Nevertheless, 
the properties of these higher- dimensional black holes would still be modified compared 
to their 4-dimensional analogues [4, 5J. 

If M* is sufficiently low, such black holes may be produced in ground-based colliders 
[6], although their appearance in cosmic rays is possible as well [7] (for reviews, see [5j[8j 
[9]). A black hole created in such trans-planckian collisions is expected to gradually lose 
its angular momentum and finally its mass through the emission of Hawking radiation 
|10j . consisting of elementary particles of a characteristic thermal spectrum, both in 
the bulk and on the brane. The emitted radiation from a higher-dimensional black hole 
created in trans-planckian collisions has been studied both analytically and numerically. 
Until recently, the Schwarzschild phase, along with a variety of additional spherically- 
symmetric black-hole backgrounds, were the most commonly studied cases. Those studies 
included the black-hole emission of lower-spin degrees of freedom [11], [T2j [T3j [TJ[ [T5j [161 
rm [T8] as well as gravitons [TH, [201 12H [22] , both on the brane and in the bulk. 

The complexity of the gravitational background around an axially-symmetric black 
hole, increased by the presence of extra dimensions, deterred many researchers from in- 
vestigating the radiation spectrum of a higher- dimensional rotating black hole. However, 
during the last two years, a plethora of studies of the Hawking radiation, emitted by 
such a black hole, appeared in the literature |23[2H[23I23I2ZII2B1I231I3Q|. These works 
offered exact numerical results that supplemented and generalized two early analytic 
works focussed on the particular case of a 5-dimensional rotating black hole [311 [32] . 
Nevertheless, to our knowledge, up to now no work has performed a complete analytic 
study of the emission of Hawking radiation from a rotating black hole in an arbitrary 
number of dimensions. 

In the present article, we consider the evaporation of a (4 + n)-dimensional rotating 
black hole into scalar degrees of freedom on the brane. We calculate the corresponding 
absorption probabilities and cross-sections obtaining analytical solutions in both the high 
and low-energy regimes. In section 2, we consider the metric corresponding to a higher- 
dimensional rotating black hole, and write down the equation for scalar fields propagat- 
ing in the projected-on-the-brane background. In section 3, we focus on the low-energy 
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regime, and solve analytically the scalar-field equation employing the matching technique 
of combining the far-field and near-horizon parts of the solution. Subsequently, we derive 
an analytic expression for the absorption probability, and produce a set of plots exhibit- 
ing its dependence on particle quantum numbers and topological properties of spacetime. 
We proceed by deriving corresponding numerical plots, and compare them to the ana- 
lytic ones, with excellent agreement. We finally derive a low-energy, simplified expression 
for the absorption cross-section, and confirm the universal behaviour characterising the 
absorption of scalar fields in this particular energy regime. In section 4, we turn our 
attention to the high-energy regime. We construct an analytic high-energy solution to 
the scalar-field equation by employing a new method involving Kummer functions, and 
calculate the absorption probability at this energy regime. By using our analytic results, 
we produce plots which we compare to the corresponding numerical ones finding again 
a very good agreement. In the same section, we determine, through numerical integra- 
tion, the exact asymptotic value of the total absorption cross-section at the high-energy 
regime. We also include an analytic treatment of the geometrical optics limit, which is 
expected to correspond to the high-energy asymptotic regime. We derive expressions for 
the absorption cross-section in this limit for three distinct kinematical cases, and compare 
them with the exact numerical results. Finally, in section 5, we state our conclusions. 



2 Gravitational Background and Field Equations 

If one accepts the prospect of the creation of a microscopic black hole during a high-energy 
particle collision, then, due to a generically non-vanishing value of the impact parameter 
between the two particles, the emergence of a rotating black hole is the most natural 
outcome. As the black hole is created in the framework of the higher-dimensional theory, 
that is characterized by a strong gravitational force, it will itself be a higher- dimensional 
object, that "feels" the extra compact, spacelike dimensions. Under the assumption that 
the black-hole horizon is significantly smaller than the size of the extra dimensions, 
the spacetime around it may be approximated by one with a single timelike dimension 
and (3 + n) non-compact, spacelike ones. A black hole living in such a background can 
have in general up to [(n + 3)/2] angular momentum parameters. However, here, we will 
be assuming that the colliding particles are restricted to propagate on an infinitely-thin 
3-brane, therefore, they will have a non-zero impact parameter only along our brane, 
and thus acquire only one non-zero angular momentum parameter about an axis in the 
brane. The background around a higher-dimensional rotating black hole with one angular 
momentum parameter is given by the following Myers-Perry solution [33] 

ds 2 = ( 1 - — dt 2 + 2a/iSin ' g dt dip - ^ dr 2 - E d6 2 

_ ( r 2 + a 2 + sm? e \ gin2 e d ^ 2 _ r% cQg2 e ^ (1) 



E r 



where 

A = r 2 + a 2 --^ T , E = r 2 + a 2 cos 2 6 , (2) 
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and df2 2 is the line-element on a unit n-sphere. The mass and angular momentum of the 
black hole are then given by 

(n + 2)A n+2 2 

M BH = — — n, J = — — M BH a , (3) 

1d7tG n + 2 

with G being the (4 + n) -dimensional Newton's constant, and A n+ 2 the area of a (n + 2)- 
dimensional unit sphere given by 

27r (n+3)/2 

A " +2 = r[(n + 3)/2] • (4) 

Since the creation of the black hole depends crucially on the value of the impact 
parameter between the two highly-energetic particles [3j, and that in turn defines the 
angular momentum of the black hole, an upper bound can be imposed on the angular 
momentum parameter a of the black hole by demanding the creation of the black hole 
itself during the collision. The maximum value of the impact parameter between the two 
particles that can lead to the creation of a black hole is j9] 



b 



max 



1 + 



n 



, (5) 



an analytic expression that is in very good agreement with the numerical results produced 
in the third paper of Ref. [3]. If we write J = bMsH/Z, for the angular momentum of 
the black hole, and use the following expression for the black-hole horizon 

that follows from the equation A(r) = 0, and the second of Eqs. ([3]), we obtain 

C = -±- • (7) 

In the above, we have defined, for convenience, the quantity a* = ajvh- Equation (JTj), 
thus, puts an upper bound to the value of the black-hole angular momentum parameter, 
which for n > 1 would have been unrestricted, contrary to the cases of n = and n — 1, 
where a maximum value of a exists that guarantees the existence of a real solution for 
the black-hole horizon. 

In this work, we will focus on the propagation of scalar fields in the gravitational 
background induced on the brane, where all ordinary particles are assumed to live. The 
4-dimensional induced background will be the projection of the higher-dimensional one 
onto the brane, and its exact expression follows by fixing the values of the additional 
azimuthal angular variables - introduced to describe the n compact extra dimensions - 
to 9i = 7r/2, for i — 2, ...,n + 1. Then, the induced-on-the-brane line-element takes the 
form 

ds 2 = (i- de + 2 T^ d dt dp - ^ 



(8) 

o f o n a 2 a sin 2 9\ 
T,d9 2 - ir 2 + a 2 + — £ — sin 2 #V- 
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Note that although the above background is very similar to the usual 4-dimensional Kerr 
one, it is not exactly the same due to its explicit dependence on the number of additional 
spacelike dimensions n. It is this dependence that will cause brane quantities to depend 
on the number of dimensions that exist transverse to the brane. 

In order to study the propagation of fields in the above background, we need to derive 
first their equations of motion. We assume that the particles couple only minimally to 
the gravitational background and have no other interactions, therefore, they satisfy the 
corresponding free equations of motion. The latter, for particles with spin 0, 1/2 and 1, 
propagating in the induced-on-the-brane gravitational background (jHJ), were derived in 
[51 [32] . For scalar fields, the field factorization 

<f>{t, r, 9, <p) = e-™ e im * R^r) Tf(0, aw) , (9) 

where T™(9,au) are the so-called spheroidal harmonics |34J, was shown to lead to the 
following set of decoupled radial and angular equations 



*( A d -^) + (^-A T )R„ em = 0, 



J_ df sine dme au) \ fm? ^ ^ e \ 
sm 9 d9 V d9 J \ sm 2 9 J 

respectively. In the above, we have defined 

K = (r 2 + a 2 )uj-am ) A™ = E™ + a 2 u 2 — 2amu . (12) 

The angular eigenvalue EJ i (aw) provides a link between the angular and radial equation. 
Its expression, in general, cannot be written in closed form, however, an analytic form 
can be found [35] in terms of a power series with respect to the parameter aw. We will 
return to this point in Section 3. 

By solving Eq. (flQl) . we determine the radial part of the field wavefunction, and, 
subsequently, the absorption probability |.4.£ )m | 2 for the propagation of a scalar field in 
the projected-on-the-brane background. This quantity appears in the differential emission 
rates for Hawking radiation emitted by the higher-dimensional black hole on the brane. 
For example, the particle flux, i.e. the number of particles emitted per unit time and 
unit frequency, has the form 

d'N _ 1 V 1 ■ |2 . , . 

dtdw ~ 2n f£ exp [k/Ta] - 1 1 lM ' 1 ' 

similar formulae may be written for the energy and angular momentum emission rates. 
In the above, k and Tjj stand for 

,_ O 171(1 T { U + l ) + { U ~ l ) a l (, A ^ 

' = "-'» n = »-^?. MTT^K ' (14) 

with Q the angular velocity of the black hole and Tjj its temperature. The absorption 
probability |^ jm | 2 depends both on particle properties (energy u, angular momentum 
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numbers £, m, etc) and gravitational background properties (number of extra dimen- 
sions n, black-hole angular momentum parameter a). As a result, it modifies the various 
emission rates from the ones for a blackbody. Equation (TlO|) has been solved analyt- 
ically only for the case of a 5-dimensional rotating black hole (and then, only in the 
low-energy regime) [32] , and numerically in [23j E2E [27] for arbitrary dimension^. In the 
next sections, we will attempt to derive analytic results, both in the low- and high-energy 
regimes, for the absorption probability for scalar fields propagating in the 4-dimensional 
spacetime of a brane embedded in the background of a rotating black hole of arbitrary 
dimensionality. 



3 Greybody Factor in the Low-Energy Regime 

In this section, we focus on the solution for the absorption probability in the low-energy 
regime. We will first derive an analytic expression by using a well-known approximate 
method. We will then plot this expression to reveal its dependence on a number of 
parameters, such as the angular momentum numbers of the particle, the dimensionality 
of spacetime, and the angular momentum of the black hole. It will also be directly 
compared with the exact numerical results derived earlier in the literature. We will 
finally derive a compact analytic expression, valid in the limit to — > 0, and comment on 
the form of the corresponding absorption cross-section and its relation to the area of the 
black-hole horizon. 



3.1 Solving the Field Equation Analytically 

In what follows, we will use an approximate method and solve first the radial equation 
of motion ( TTOT) at the two asymptotic regimes: close to the black-hole horizon (r ~ r\), 
and far away from it (r ^> r/J. The two solutions will then be stretched and matched 
in an intermediate zone to create a smooth analytical solution extending over the whole 
radial regime. 

We first focus on the near-horizon regime and make the following change of variable 



= 2 i 2 T = ^ ~ f) T 2 i 2 ' ( 15 ) 

r z + cr ar r A + cr 



where, for convenience, we have defined the function Air) = (n + 1) + (n — 1) a 2 /r 2 . 
Then, Eq. ( FlUj) . near the horizon (r ~ r^), takes the form 



d 2 P dP 
/(l- /) ^ + (l-^/ ) - + 



Kl AT (1 + al) 



Alf(l-f) A*(l-f) 



P = 0, (16) 



where now 

A* = (n + 1) + (n — 1) al, K* = (1 + al) lu* - a*m , (17) 



lr The angular equation pip was also solved numerically in |25j . and the angular distribution of the 
emitted radiation was found. 
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with = ur h . We have also defined the quantity 

n(l + a 2 ) Aal , , 

d. = i+ 4: if- (18) 

By making the field redefinition P(/) = / a (l — f) F(f), Eq. ([TBI) takes the form of 
a hypergeometric equation [36J: 

r/ 2 F rfF 

/(W)^ + [c-(l + a + &)/] — -a&F = 0, ( 19 ) 

with 

a = a + P + D*-l, b = a + j3 , c=l + 2a. (20) 

The power coefficients a and (3 will be determined by two constraints (following from the 
demand that the coefficient of F(f) is indeed — ab) that have the form of second-order 
algebraic equations, namely 

« 2 + ^ = 0, (21) 

and 

/3 2 + m-2) + ^- ^ { \t =°- ( 22 ) 

The general solution of the hypergeometric equation ( 1T91) . combined with the relation 
between P(f) and F(f), leads to the following expression for the radial function P(f) 
in the near-horizon regime: 

P NH (f) = A„r(l-ffF(a,b,c; f) 

+ A + f- a (l- ffF{a-c+ 1,6 -c+ 1,2- c; /) , (23) 

where A±. are arbitrary constants. Solving Eq. (I2~TI) . we obtain the solutions 

a ± = . (24) 

Near the horizon, r — > and /(r) — > 0. Then, the near- horizon solution (1231) reduces to 

PnhU) ^ A- / ±iK */ A * + A + f^ iK *l A * = A. e ±iky + A + e Tlky , (25) 

where, in the second part, we have used the definition for k given in Eq. (fTl"]) . and the 
tortoise-like coordinate y = r h {l + a 2 ) ln(/)/A*. Note, that although the coordinate y 
is not identical to the usual tortoise-one, defined by dr^/dr = (r 2 + a 2 )/A(r), used in 
Kerr-like backgrounds [21], it holds that 

dy _ f A\ (r 2 + q 2 ) 2 /r h y-z dr, 

dr \Aj (r 2 + a 2 ) 2 \ r J dr ' { ' 

Therefore, in the limit r — > r^, the two become identical, and the near-horizon asymptotic 
solution assumes, as expected, the free-wave form in terms of the tortoise-coordinate 
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Imposing the boundary condition that no outgoing mode exists near the 
horizon, we are forced to set either A_ = or A + = 0, depending on the choice for a. 
As the two are clearly equivalent, we choose a = <x_, and set A + = 0. This brings our 
near-horizon solution to the final form 

P NH (f) = A.f a (1 - ff F(a, b, c; f) . (27) 
We finally turn to Eq. (|22|) for the (3 power coefficient. This admits the solutions 

1 



* 2 



(2f 



The sign appearing in front of the square root will be decided by the criterion for the 
convergence of the hypergeometric function F(a, b, c; /), i.e. Re (c — a — b) > 0, which 
demands that we choose (3 = 

We now turn our attention to the far-field regime. Making the assumption that 
r ^> Th, and keeping only the dominant terms in the expansion in terms of 1/r, the 
radial equation (Tl0|) takes the form 

d 2 P 2dP ( 9 Er + a 2 uo 2 \ 

+ r— + K - — J p (r) = • ( 29 ) 



dr 2 r dr 



where Ef 1 is the angular eigenvalue. The substitution P{r) = P(r) brings the above 

into the form of a Bessel equation for P(r), and the overall solution in the far-field limit 
can be written as 

P FF (r) = ^=J» (wr) + ^pY v (ur) , (30) 

where J u and Y v are the Bessel functions of the first and second kind, respectively, with 
v = a/ P\ a + a 2 uj 2 + 1/4, and B\ 2 integration constants. 

In order to construct an analytic solution extending over the whole radial regime, 
we need to smoothly match the two asymptotic solutions, derived above, in some inter- 
mediate regime. Before doing so, we first need to extrapolate (stretch) the two solutions 
towards this regime. To this end, we shift the argument of the hypergeometric function 
of the near-horizon solution from / to 1 — / by using the relation 



p NH (f) = A_r(w) 



r(c)T(c-a-b) , , 

— — F(a, b, a + b - c + 1; 1 - /) 

r(c — a) T(c — b) 



T(c)T(a + b-c) 



+ (l-/r a - 6 w p ;. m . L F (c-a,c-b,c-a-b+l;l-f) 

r(a)r(6) 

The function f(r) may be alternatively written as 

A* 1 f r h\ + 



.(31) 



/(r) = l--f T ^— ^ = 1- - \ ' f 2 , (32) 

r n 1 r ^ _|_ V r / \ r l r h) + a * 
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where we have used the horizon equation, A(r^) = 0, in order to eliminate fi from the 
above relation. In the limit r ^> r^, the (r/r/J 2 in the denominator of the second term 
above is dominant, and the whole expression goes to unity for n > 0. Thus, in the limit 
/ — > 1, the near-horizon solution (13T|) takes the form 

P NH (r) ~ A ir -( n+1 ^ + A 2 r^^ +D ^ , (33) 



with 



- A_[(l + ajr, j r(c _ fl)r(c _ 6) 

4 4 m i „2 W n+i r (/?+fl.-2) r( c )r(a + &-c) 

^ - A.L(l + a,)r fc j r . (34) 



Next we need to stretch the far-field asymptotic solution (I30I) towards small values 

Then, in the li 

Bi ^ j v /jB r+ a2 ^ 2 + i /4 



of the radial coordinate. Then, in the limit oor — > 0, we find 



P F F{r 



y/rT (y£™ + a 2 cj 2 + 1/4 + 1 

" ^ r v v Er + aw +1/4 J \y) 



• (35) 



We notice that both 'stretched' asymptotic solutions have reduced to power-law 
expressions in terms of the radial coordinate r, however, the different power coefficients 
prevent the exact matching. In order to overcome this obstacle, we will expand these 
power coefficients in the limits (cur^) 2 <C 1 and (a/r^) 2 <C 1. It is these approximations 
that will limit the validity of our result for the absorption probability to the low-energy 
and low-angular-momentum regime. Note, however, that in order to improve the accuracy 
of our result, no approximation will be made in the arguments of the Gamma functions 
involved in Eqs. (1331) and (|35|) . In order to follow the aforementioned line of action, we 
need to know the analytic expression for the angular eigenvalue E™. According to [35], 
this may be written as a power series in terms of the parameter (ouj), namely 

oo 

K = E ^ (*")" • (36) 
rc=0 

For the purposes of our analysis, we have calculated the coefficients f l ™ up to fifth order, 
and found the result 

,m , ,s , , x 2 [2™ 2 - 2/ (/ + !) + !] 



E™ = /(/ + !) + (aw) 



+ (aw) 4 



+ 



(21 - 1) (21 + 3) 
2 [-3 + 17/ (/ + !) + l 2 (l + 1) 2 (2/ - 3) (21 + 5)] 



(21 - 3) (21 + 5) (21 + 3) 3 (2/ 



31(1 + V 



(21- l) 2 (2/ + 3) 2 [(21- 1) (2/ + 3) (2/ -3) (2/ + 5) 



2m 4 [48 + 5(2/-!) (21 + 3)] , 
(2/ -3) (2/ + 5) (2/ -1)3(2/ + 3)3 l "" 



S 



with f[ m = /| m = = 0. The above form will be used at every place where E™ appears 
in Eqs. (1331) and (1331) . The only exception will be in the power coefficients, where terms 
of order (au) 2 , or higher, will be ignored. Following these assumptions, the two power 
coefficients in Eq. (|33|) reduce to 



-(n+l)P ~ l + 0(Lulaia*LU*) , 
(n + + A, -2) ~ -(Z + l) + 0(a&c&a^). 
while the power coefficient in Eq. (1351) takes the form 



E ra + G 2^2 + 1/4 ^ ( / + _ j + . 



(38) 



(39) 



By using the above results, one can easily show that both Eqs. ( 133]) and (I35I) reduce 
to power-law expressions with the same power coefficients, r l and r~" +1 >. By matching 
the corresponding coefficients, we determine the integration constants B\ 2 in terms of 
the other parameters of the theory. Their ratio, which, as we shall see, appears in the 
expression of the absorption probability |^4i, m | 2 , is found to be 



B = 



Bi 
Bo 



21+1 



E™ + a?u 2 + 1/4 



x 



* \ur h (l + a2)n+i 

T 2 L/Ef + a 2 u 2 + 1/4) r(a + P + D, - 1) r(a + (3) T(2 - 2/3 - D») 
r(2/3 + D, - 2) T(2 + a - (3 - A) T(l + a - /?) 



.(40) 



The last step in our calculation of the absorption probability involves the expansion 
of the far-field solution (!30j) to infinity (r — > 00). This leads to 



P FF {r) 



V2 



[Bx + iB 2 ) 



+ 



(Bi - iB 2 ) 



A 



(00) 



1(00) ^_ 



+ A 



(41) 



As expected, at large distances from the black hole, where the effect of the angular 
momentum parameter a is almost negligible, the solution for the scalar field reduces to a 
spherical wave, as in the Schwarzschild case [TU [231 [25] . Then, the absorption probability 
is defined through the expression 

2 



Lm I 



1 



.4 



(00) 

out 



A 

B 



(00) 



1 



B 1 - iB 2 



B 



B x + iB 2 
2i (B* - B) 



BB* + i(B*-B) + l' 



(42) 



The above result, together with the expression for B given in Eq. (|40p . is our main 
analytic result for the absorption probability for scalar fields valid in the low-energy 
and low-angular-momentum regime. It can be easily checked that it reduces smoothly 
to the corresponding result for a scalar field propagating in a Schwarzschild-like, higher- 
dimensional background projected onto the brane [TTJ, if we set a = 0. 
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Figure 1: Absorption probability |^4; jm | 2 for brane scalar particles, for n = 2, a* = 0.4, 
and for the modes (from left to right) (7 = m = 0), (I — 1, m — —1,0, 1), and (/ = 2, 
m = —2, —1,0, 1,2). The solid lines correspond to our analytic results, and the dashed 
lines to the exact numerical ones. 

3.2 A Comparison with the Exact Numerical Solution 

In this section, we proceed to study in detail the properties of the absorption probability 
|*4;, m | 2 , as this was derived above by using purely analytic arguments. To this end, we will 
plot our main result, given by Eqs. fj4"2l and fj4*0|) . as a function of the energy parameter 
coT/j and for a variety of values of the other parameters of the theory, namely the angular 
momentum numbers (l,m) of the scalar particle, and the topological parameters (a*,n) 
of the spacetime. At the same time, a direct comparison of our low-energy analytic 
result to the exact numerical value for |^4; jm | 2 - derived in [91 [23j [25] for the purpose 
of computing the Hawking radiation emission spectra - will be performed in order to 
examine the range of validity of our approximations. 

In Fig. [1], we plot the absorption probability |*4; im | 2 for scalar particles, for fixed an- 
gular momentum of the black hole (a* = 0.4) and number of extra dimensions (n — 2), 
and for a variety of modes with different angular momentum numbers (l,m). Through- 
out our paper, our analytic results will be plotted by using solid lines while the exact 
numerical results, reproduced from [9J [231 125] , will be denoted by dashed lines. In Fig. [lj 
both sets of lines are shown for all modes, and the agreement between them at low energy 
is indeed remarkable. Although a small deviation appears when the energy parameter is 
taken beyond the low-energy regime, the qualitative agreement between the two sets of 
results remains excellent. 

Focusing now on the dependence of the absorption probability on particle param- 
eters, we observe that, similarly to the Schwarzschild case [TT] . it is again the lowest 
partial wave, with / = 0, that dominates in the low-energy regime, with all higher modes 
increasingly suppressed. This behaviour is valid for all values of a* and n, as long as 
attention is focused on the low-energy regime. Looking next at the relative behaviour of 
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Figure 2: Absorption probability |*Ai, m | 2 for brane scalar particles, (a) for the modes 
(I = l, m = -1,1), for n = 2 and a, = (0,0.2,0.4,0.6,0.8), and (b) for the modes 
(I = l,m = —1,0,1), for a* = 0.4 and n = (1,2,4,7). Again, the solid lines correspond 
to our analytic results, and the dashed lines to the exact numerical ones. 



modes with the same angular momentum number I but different number m, we easily 
note that the modes with m < are the dominant ones with increasing m causing sup- 
pression. We also observe that, for modes with m < 0, the absorption probability always 
remains positive, while, for modes with m > 0, a negative- valued region for |^ jm | 2 always 
appears in the low-energy regime. The latter effect is due to the so-called superradiance 
[38] - the enhancement of the amplitude of an incoming wave by a rotating black hole 
that results in a reflection probability larger than unity, and thus a negative absorption 
probability according to |^4; jm | 2 = 1 — |7£/ im | 2 . We will say a little more on the origin of 
this effect in the following subsection. 

We now turn to the dependence of the absorption probability on the parameters of 
spacetime. In Fig. [21 we depict the behaviour of |*4.z, m | 2 in terms of the angular momentum 
parameter a* and number of extra dimensions n. We have chosen the indicative case of 
/ = 1, and, in Fig. [2](a), depict the behaviour of the m = — 1 (from central to left) and 
m = 1 (from central to right) modes, for a range of values of a*. Then, in Fig. [2]^b), 
we present all three modes, with m = 0,±1, for various values of n. As before, both 
the analytic, low-energy results as well as the exact numerical ones are shown. Again, 
the agreement between the two sets of results in the low-energy regime is remarkably 
good. The alert reader may note that, in general, the agreement between the two sets is 
improving as n increases: this is due to the fact that, according to Eq. (132]) . an increase 
in the number of extra dimensions improves the accuracy of the assumed behaviour of 
the function f(r) at infinity, and consequently our approximation. In addition, terms 
that have been neglected during the matching of the two asymptotic solutions under the 
low-energy and low-a* assumption, such as the K^/Al in Eq. (128]) . become even smaller 
for large values of n, thus improving the accuracy of our analysis. 

According to Fig. [2](a), for fixed n, the non-sup err adiant modes, with m < (al- 
though not shown, the m = mode exhibits a behaviour similar to the m = —1), are 
enhanced with the angular momentum of the black hole in the low-energy regime, while 
the superradiant modes, with m > 0, are suppressed both in the superradiant and non- 
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superradiant energy regime. Turning to Fig. |2](b), we note that, for fixed a*, an increase 
in the number of extra dimensions n also leads to an enhancement of the value of the 
absorption probability for the non-superradiant modes. For the superradiant ones, the 
behaviour of |*4.z, m | 2 depends on the energy regime we are looking at: while it is sup- 
pressed in the superradiant regime, it is enhanced in the non-superradiant one. From 
both figures, it becomes obvious that the superradiant effect becomes more important as 
either a* or n increases. This enhancement, for brane-localised scalar particles, both in 
terms of the angular momentum of the black hole and the dimensionality of spacetime 
was also noted in the literature [231 I2U [28] . 



3.3 The Low-Energy Asymptotic Limit of the Cross-Section 



Our last task regarding the behaviour of the absorption probability |^4i iJn | 2 in the low- 
energy regime will be the derivation, from Eqs. (l4"0j) and fj42l . of a compact analytic 
expression valid in the limit u — > 0. This simplified analytic expression will be used 
to explain some of the features discussed in the previous subsection. In addition, from 
this, the asymptotic low-energy value of the corresponding absorption cross-section for 
scalar fields in the background of a projected-on-the-brane rotating black hole will also 
be determined. 



We start our analysis by noticing that, according to Eq. (|4"CT|) . in the limit u) 
B ~ u- {2l+ V and, therefore, BB* » i(B* - B) > 1. Then, Eq. simplifies to 



2 _ 2i(B* - B) _ 



1 

~B* 



BB* \B 

Substituting for B using Eq. ( )40i) . and the fact that a is purely imaginary, yields 



0. 



(43) 



l.m I 



-2in (ujr h /2) 



21+1 



T(2(3 + D, - 2) 



(/ + !)!*(/ + I) (1 + 



2Z + 1 

iiy—T{2-2(3-D*) 



x 



1 



X 



T(a + (5 + D* - 1) F(-a + (3 + D # - 1) T(a + (3) T(-a + (3) 
F(2 + a-(3- A) r(-a + p + D* - 1) T(l + a - (3) T{-a + (3)- 

r(2 - a - (3 - D # ) T{a + (3 + D,~ 1) r(l - a - (3) T{a + (3) 
Si x S 2 x S 3 , 



(44) 



where Si, S 2 and S3 are defined by the quantities on each of the three lines above. 
Focusing first our attention to S3, and using the Gamma function relation T(z)T(l — z) = 
7t/ sin7rz [36], this can be written as 

-7T 2 sin(27ra) sm%(2P + D*) 

3 sin ir(a + (3 + D*) sin 7r(— a + (3 + D*) sin ir(a + (3) sin 7r(— a + (3) 

From the factor u 2l+1 in Eq. (j4"4"j) . it becomes clear that the expression for the absorption 
probability at the very low-energy regime is dominated by the lowest partial waves, a 
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property that is in accordance with the results presented in the previous subsection. 
Then, assuming that m is small and a* < 1, the limit u — > is equivalent to a — > 0. 
Expanding terms in S3 and £2 to linear order in a gives 

2n 3 a wnirpp + D„) _ 1 

(46) 



sin 2 n{(3 + D,) sin 2 vr/T T(J3 + D* - l) 2 T((3) 2 ' 

Using the additional Gamma function relation T(z)T(—z) = — Ti/zsimrz allows the over- 
all expression for the low-energy limit of the absorption probability to be written as 

|2 _ 4tt (ur h /2) 2l+1 K* sin 2 n(2(3 + Dj T 2 (2/3 + g, - 2) T 2 (l - g) (2 - g» - 2(3) 

7Th\ — 2Z4-1 * 

A, (1 + alY^r (I + I) r 2 (/ + I) r 2 (/3 + A, - 1) sin 2 tt(/3 + £>*) 

(47) 

In the above, we have also used the definition a = —iK*/A#, where, from Eq. ( fl~7l) . 

X* = (1 + a 2 ) - a*m = r h (l + a 2 ) (w-mfi). (48) 

By using Eq. (f2"8"j) . one may easily conclude that the quantity (2 — — 2/?) is always 
positive, while the same also holds for A* = (n + 1) + (n — l)a 2 , for all values for 
a* and n > 0. Therefore, the overall sign of |^ jm | 2 is determined by the sign of K*, or 
equivalently of (u—m fl). A negative sign for the latter combination arises only for m > 0, 
and denotes the occurrence of superradiance, with |^4; jm | 2 acquiring a negative value. 
The superradiance domain arises in the low-energy regime and extends over the range 
of values < uo < uo s = m Q. The larger the value of the angular momentum parameter 
a*, the larger the rotation velocity Q of the black hole, and thus the more extended the 
superradiance regime becomes. This is indeed in agreement with the behaviour found in 
the previous subsection. 

In what follows, we focus on the dominant s-wave with I = m = 0. As we will see, 
this will be the only partial wave with a non-vanishing low-energy asymptotic value of 
the absorption cross-section. In order to simplify further Eq. ( I47p . we need also to expand 
the expression of (3 in the limit lu — > 0. It is easy to see that, in this limit, (3 = + 0(lu 2 ), 
which then allows us to write 

1 01 A* (1 + a 2 )- 1 /^!) (2 - D m ) { ' 

The corresponding absorption cross-section for the dominant partial wave is then given 

by m 

°o = - \A \ 2 = 4vr (r 2 + a 2 ) [ J *> + ... . (50) 

uj 2 [(n + 1) + (n - 1) ai\ ( 2 _ D *) 

According to the above result, the absorption cross-section for the lowest mode / = 
reduces to a non-vanishing asymptotic value, as uj — > 0. We may easily see, from Eq. 
(j47p . that the low-energy behaviour of the absorption probability for any higher partial 
mode will be governed by the factor uj 21+2 , thus leading to an absorption cross-section 
proportional to u 21 . Therefore, for all partial waves with 1^0, the partial cross-section 
goes to zero, as uj — > 0. 
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Figure 3: Absorption cross-section o"o (in units of r^) for the lowest scalar mode I = 0, 
for n = 2 and a* ranging between and 0.6. As before, the solid lines correspond to our 
analytic results, and the dashed lines to the exact numerical ones. 



In the case of scalar particles propagating in a Schwarzschild-like projected-on-the- 
brane line-element, the low-energy asymptotic value of the absorption cross-section of 
the lowest, dominant partial wave was shown to be equal to the horizon area of the 4- 
dimensional black hole, 47rr 2 , regardless of the number of extra dimensions pJ],[T3]. We 
would like to demonstrate that a similar relation holds in the case of an axially-symmetric 
brane background. According to Eq. fl50|) . <7o is indeed proportional to the horizon area 
of the four-dimensional rotating black hole ^{r\ + a 2 ), however, the relation involves a 
multiplicative factor which is both (a*, n)-dependent. In Fig. [31 we plot cr , for fixed n 
and various values of a*, by using both our analytic expression and the exact numerical 
result. The latter set of results reveal that the asymptotic low-energy cross-section is 
indeed equal to the horizon area of the black hole, regardless of the number of extra 
dimensions. For small values of a*, the multiplicative factor appearing in Eq. (|50|) is very 
close to unity and our analytic expression closely reproduces the exact numerical one. 
As a* increases though, the range of validity of our approximation is exceeded and a 
deviation starts appearing, as expected. If we keep a* fixed and vary n instead, a similar 
behaviour appears, with values of n < 1 leading, in general, to a value smaller than the 
exact result, and values of n > 2 to a value larger than the exact result. The magnitude 
of the deviation depends again on the value of a*. The above results demonstrate the 
universal behaviour for the lowest partial mode of a scalar particle according to which its 
partial cross-section equals the area of the black-hole horizon in the low-energy regime. 
This result holds not only for a spherically-symmetric brane line-element but also for an 
axially-symmetric one, and is independent of the number of the transverse-to-the-brane 
spacelike dimensions. The latter result was also reproduced for the particular case of a 
5-dimensional bulk in [30J. 
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4 Greybody Factor in the High-Energy Regime 



In this section, we turn our attention to the high-energy regime, and present a way to 
compute the absorption probability in this region by using again the matching tech- 
nique described in section 3. For this, we are going to use the near-horizon solution we 
have already computed in Section 3.1, and we shall construct an approximate, high- 
energy, far-field solution that will allow us to do the matching solely in the high-energy 
regime. Next, we are going to compare the absorption coefficient produced this way with 
the exact numerical results of the literature. Finally, generalizing well-known results for 
Schwarzschild black holes, the geometrical optics limit value of the absorption cross- 
section will be computed, and its connection to the high-energy asymptotic value found 
by numerical analysis will be investigated. 



4.1 Analytic Construction of the Solution 

In order to construct a solution over the whole radial domain that will be valid in the 
high-energy regime, it is necessary to do the matching of the stretched near-horizon 
and far-field asymptotic solutions without resorting to the low-energy approximation 
{urh <C 1) employed in section 3.1. To this end, we will try to find a new far-field 
asymptotic solution that will satisfy the field equation only in the high-energy limit. 
Moreover, the exact form of this solution should be such that, when stretched towards 
small values of the radial coordinate, it reduces to a power-law expression with identical 
power coefficients to the ones appearing in the stretched near-horizon solution, thus 
allowing for a perfect matching in the intermediate regime. We remind the reader that 
the stretched near- horizon solution, (Eq. I331) . was found to be of the form 

P NH ~ A ir -( n+ W + A 2 r( n+1 W +D *-V = A ± r^ + A 2 r^ . (51) 

The differential equation that our far-field solution needs to satisfy is Eq. (TTOl) in the 
limit r ^> rh, or equivalently 

^' + ^,^ + d-S±4U, W = 0, (52) 



dz 2 z 2 + a\ dz \ z 2 + a\ 

where we have made the change of variable z = ur, and defined for convenience a\ = au. 
Let us consider the following trial, special solution to the above equation 

P FF = e'^'r^M (1 +/3i,2 + 2ft; 2iur) , (53) 

where M(a,b;w) is the first Kummer function [36]. Solving Eq. fl53|) for M(d,b;w) and 
substituting in the confluent hypergeometric equation, that the Kummer functions sat- 
isfy 

d 2 M , t dM , . 

_ + ( l_„ ) __ aMW = 0, (54) 
with a = 1 + ft, b — 2 + 2ft and w = 2iur, we finally obtain the equation 

d 2 P FF 2 dPpp ( ft (ft + m . . 

-Jz^ + z — + { 1 ~ 2 PfF = ' (55) 
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Subtracting Eq. (151)1) from Eq. ([52]) . we find 



From the above, we may conclude that, for large z = cur, the constructed solution ( 1531) 
satisfies the far-field (r ^> r^) equation (I52I) up to second order in (1/z). Therefore, in 
the limit u — * oo, our trial solution (|53|) is indeed a good approximation to the exact 
solution of the field equation in the far-field regime. 

By following the same method, a second special solution to Eq. (152]) may be con- 
structed that has the form 

P FF = e-'^'r'^- 1 U {-fa, -2/32; 2za;r) , (56) 

where £7(a, 6; u>) is the second Kummer function. The second special solution follows 
from the first one under the replacement fa — > —1— fa, which preserves the structure of 
Eq. ( 1551) . Therefore, this solution, too, is a good approximation to the exact solution of 
Eq. (152|) in the limit ui — > oo. Therefore, the total solution in the far-field domain, valid 
only in the high-energy regime, takes the form 

P FF = B + e~ iu>r r Px M{fa + 1, 2fa + 2, 2iur) + 5_ e~ iuir r " /3 ^ 1 U{-fa, -2fa, 2iur) . (57) 

The Kummer functions have been used before in the construction of the far-field solution 
of a general spin-s field [21 [11], however, the projected-on-the-brane background in that 
case was that of a spherically-symmetric Schwarzschild-like one. We may easily check 
that, in the limit a — > 0, the above solution reduces to the one for a scalar field propa- 
gating in the far-field domain of a spherically-symmetric brane black-hole background. 

As an additional independent check of the above analysis, one may observe that, 
under the demand that Eq. (155!) matches exactly Eq. (|52l) . we are led to 

fa = ~ + y/E? + a 2 LJ 2 + 1/4 = -- + u . (58) 

But then, since a = v + 1/2 and b = 2v + 1, the first Kummer function reduces to 
M{y + 1/2, 2v + 1; 2iur) ~ e lujr r~ v J v {ur) [36J. By using this result - as well as a similar 
one for the second Kummer function U - in conjunction with the relations ( |53il and 
fl56l) . we may see that the far-field solution (1301) . found in section 3.1, is duly restored as 
expected. 

Nevertheless, in this section, the approximate solution (l57|l will be used instead, since, 
as was mentioned earlier, its use will allow us to achieve a perfect matching between the 
stretched near- horizon and far-field solutions - a feature that, as we saw, was not possible 
when the asymptotic solution (1301) was used instead. To this end, we stretch Eq. (1571) to 
small values of ur [36J to obtain 

P FF = B+r f3i +B r P2 r( ~^ 2 ~ X) (2zu) 2 ^ +1 . (59) 

r(-p 2 ) 
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As we hoped, the stretched far-field solution contains powers of r that exactly match the 
ones appearing in the stretched near- horizon solution (lo"Tj) . Then, by matching also the 
corresponding multiplicative coefficients, we find 



B 



B_ 



[(1 + <£) r n h +1 ] 



n+11 2-2/3-A. 



X 



(2iu) 2 ^+ 1 

T{2/3 + D*-2)T(2 + a-/3- D*) Y{1 + a - (3) T{-(3 2 ) 
T(a + P + D*- 1) T(a + (3) T(2 -2(3- D*) T{-2(3 2 - 1] 



(60) 



We remind the reader that the coefficients -D*, a and (3 are given in Eqs. ( |T8i) . ( 12^1) and 
( |28i) . respectively, while the coefficients (3\^ are defined in Eq. ( I5TI) . 

In order to finally compute the absorption probability, we need first to expand the 
far-field solution (1571) in the limit r — > oo. Then, we find 136 



r(2 + 2/3 x ) e 



«7r(/3i+l) 



r(i + /3i) (2^)^ +l 



,02 



+ 



s+ r(2 + 2A 



r r(l + ft) (2io;)ft 



+i 



A- 



(oo) e" 



+ A (°°) e _ 

' -"-out 



+ ... 



(61) 



As we see, the far-field solution ( 157|) also reduces to a spherical free-wave solution in 
the asymptotic infinity. We may thus use once again the standard definition for the 
absorption probability to determine its value 



IA 



l - 



A, 



(oo) 
out 



A 



(oo) 



r(2 + 2p 1 



r(2 + 2ft) e i -(ft+ 1 ) + £ r(/3i + 1) (2iu)K 



+/32+1 



(62) 



The above expression, combined with Eq. (1601) . gives the absorption probability for scalar 
fields, valid only in the high-energy regime, but with no restrictions on the value of the 
angular momentum parameter a apart from the upper bound of Eq. (J7|). The correspond- 
ing absorption cross-section, valid in the high-energy regime, then follows by using the 
formula a^ m = ir\Ai >m \ 2 / u 2 and Eq. (I6"2l . 

In Figs. H](a,b), we depict the absorption probability |A jm | 2 and cross-section a^ m , 
respectively, for brane scalar particles, in the high-energy regime. Once again, the solid 
lines correspond to our analytic results, following from Eq. ( 1621) . and the dashed lines to 
the exact numerical ones. Figure H](a) depicts the absorption probability for the indicative 
case of the three lowest partial modes. We notice that our analytic results match the 
exact numerical ones for large enough value of the energy to. The lower the value / of the 
partial mode, the sooner the two results coincide. As the energy parameter increases, the 
absorption probability quickly tends to unity - as expected, highly energetic particles 
always overcome the gravitational barrier outside the horizon of the black hole. This 
asymptotic behaviour, although successfully reproduced by our analytic results, does not 
allow us to appreciate the agreement between the two sets of results. For this reason, in 
Fig. IHb), we plot the absorption cross-section for the three lowest partial waves summed 
over m: in the high-energy regime, a^ m ~ 1/uj 2 for all modes, and the asymptotic regime 
is significantly more extended. The agreement between our analytic and numerical results 
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Figure 4: (a) Absorption probability \Ai\ 2 for brane scalar particles, for the modes / = 
0, 1, 2 and m = 0, for n = 2 and a* = 0.3, and (b) absorption cross-section 07 (in units 
of r|), for the modes I = 0, 1, 2 (summed over m), and for the same values of n and a* 
as before. Again, the solid lines correspond to our analytic results, and the dashed lines 
to the exact numerical ones. 



is now much clearer: for the mode I = 0, the two results completely coincide, while, 
for I = 1,2, the exact matching is achieved at gradually larger values of energy. Once 
matched, the two results remain identical as u increases further until the zero asymptotic 
value - for the individual partial modes - is reached. Given the increased difficulty in 
integrating numerically both the radial and angular part of the scalar equation of motion 
over an extended energy regime (for example, see [T3J [22] )> the above solution could be 
used to analytically extrapolate a numerical solution to arbitrarily large values of the 
energy parameter uj. 

4.2 High-Energy Asymptotic and Geometrical Optics Limits 

As the energy of the particle emitted from a black hole increases, the total absorption 
cross-section o" a b s = J2i m a i,m reaches a high-energy asymptotic value, in an oscillatory 
way. Although each partial cross-section a^ m asymptotes zero at the high-energy regime, 
the superposition of an infinite number of partial waves, each one reaching its maximum 
value at a gradually larger value of u as / increases, creates this constant asymptotic 
value. This asymptotic limit has been studied in the past for a Schwarzschild black hole, 
both in the four-dimensional [4UI |4T| |4"2"| 143"] and (4 + n)-dimensional case (T3J [44]. For a 
rotating black hole, the corresponding study was performed in 4 dimensions in |45j . and 
in 5 dimensions in [30J. 

Here, we will attempt to give a comprehensive study of the high-energy asymptotic 
limit of the total absorption cross-section for scalar fields living on the brane-induced 
line-element of a (4 + n)-dimensional rotating black hole. As we will see, similarly to the 
case of a Schwarzschild-like induced-on-the-brane line-element, the number of transverse 
dimensions, although inaccessible to the brane-localised scalar fields, affects the value 
of the high-energy asymptotic limit of the absorption cross-section. The value of the 
angular momentum parameter a* of the higher- dimensional black hole will also be found 
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Figure 5: Partial absorption cross-sections 07, for the modes I = 0, 1,2,3,4 (lower set of 
curves), and the total absorption cross-section cr a b s (upper curve) for n = 2 and a* = 0.3, 
both in units of r\. The dashed line denotes the value obtained by using the geometrical 
optics limit fl69|) . 



to have an effect on the value of cr a b s . Although our analytic results describe successfully, 
as we saw in previous sections, both the low- and high-energy regimes of the absorption 
probability and cross-section, no analytic solution currently exists that smoothly connects 
the two solutions over the intermediate energy regime. The emergence of the high-energy 
asymptotic limit of the total absorption cross-section strongly relies on the contribution 
of the low-energy regime (where the I = mode dominates), the intermediate-energy 
regime (where all modes have a significant contribution) and the high-energy regime 
(where higher modes dominate). As a result, in order to accurately derive the high- 
energy asymptotic value of cr a b s , and in the absence of a global analytic solution, the use 
of exact numerical analysis is imperative. 

In Fig. we therefore present exact numerical results for the absorption cross-section 
o i = X]m a ': m f° r ^ ne partial modes I = 0,1,2,3,4, as well as for the total absorption 
cross-section cr a b s for a brane-localised scalar field, for the indicative case of n = 2 and 
a* = 0.3. The emergence of a constant high-energy asymptotic value for the total cross- 
section is obvious. In the context of our analysis, we have studied the behaviour of o" a b s 
for a range of values of n and a*, with our results being displayed in Table 1. Note that, 
while for low values of n and a*, a relatively small number of partial modes needs to 
be summed over (for the case depicted in Fig. [5l only modes up to / = 5 were included 
in the calculation), as either a* or n takes large values, an increasing number of partial 
waves needs to be taken into account. In addition, as either a* or n increases further, the 
asymptotic value of cr a b s emerges at continuously larger values of lu, which significantly 
increases the computing time. 

From the entries of Table 1, one may observe the strong dependence of the high- 
energy asymptotic limit of the absorption cross-section on both the number of transverse- 
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a„\n 



1 



2 



3 



4 



5 



6 



0.0 
0.3 
0.5 
0.7 
1.0 
1.5 
2.0 



12.6 
13.6 
15.7 
18.7 
25.1 
40.7 
62.8 



9.6 

10.4 

11.5 

13.2 

16.6 

24.1 

33.6 



8.2 

8.6 

9.5 

10.7 

13.2 

18.6 

25.2 



7.3 

7.6 

8.4 

9.4 

11.4 

15.8 

21.1 



6.7 

7.0 

7.6 

8.5 

10.3 

14.0 

18.7 



6.2 

6.5 

7.1 

7.9 

9.5 

12.9 

17.2 



Table 1: High-energy asymptotic values of the total absorption cross-section cr a bs, in units 



to-the-brane dimensions and the angular momentum of the black hole. As in the non- 
rotating case [13], cr a b s is strongly suppressed as n increases. On the other hand, an in- 
crease in the value of a* causes an enhancement in the value of cr a b s - For a* = 0, the values 
of cr a bs match, as expected, the ones obtained for a scalar particle in a Schwarzschild-like 
projected brane background [13]. We would also like to note that a feature that seemed 
to hold in the 5-dimensional case [30], namely that the high-energy asymptotic value of 
cr a bs is close to the low-energy one, disappears for general n. This can be clearly seen in 
Fig. [5j or by comparing the entries of Table 1 with the value of do, from Eq. (1501) . that 
dominates the low-energy value of a a b s : for general n, the two sets of values are distinctly 
different. 

In the case of a non-rotating black hole, the geometrical optics limit has been success- 
fully used to explain the high-energy asymptotic value of the absorption cross-section 
dabs both in the pure 4-dimensional case [HH HTJ H2J H3] and (4 + n)-dimensional one 
[T3JI33]. In the higher-dimensional case, for particles living on the brane, the geometrical 
optics analysis showed that the Schwarzschild black hole behaves as a perfect absorber 
of a radius given by 



The absorption cross-section is then given by the target area, <T a bs = ^ r 1- The values 
following from this expression, for different n, are in perfect agreement with the numerical 
ones found in [13] , and displayed here in the first row of Table 1 . 

Here, we will attempt to perform a similar study, in an axially-symmetric black- 
hole brane background, and investigate the potential connection between the analytic 
values that follow from this analysis and the exact numerical ones depicted in Table 1, 
for a* ^ 0. For this, we will closely follow the method described in [45J. Although that 
formalism was developed for the case of a pure 4-dimensional Kerr black hole, it holds 
identically for the case of a projected-on-the-brane rotating black hole, with the only 
difference appearing in the exact expression of the metric function A(r). We therefore 
present here only the basic assumptions and the final equation that describe the particle's 
trajectory. 



of r^, as a function of n and a*. 




(63) 
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The line-element (jSJ), in which the brane-localised particles propagate, is invariant 
under translations of the form t —>■ t + At and <fi — > <ft + A<f>. The corresponding Killing 
vectors £^ = 5f and \ = then lead to the conserved conjugate momenta p t = —E 
and = — Z^. The brane metric also possesses a Killing tensor that leads to an 
additional conserved quantity Q = ^uP^p u — (Ea + L z ) 2 . Combining the above, the 
equation of motion of a particle with rest mass m, i.e. p^p^ = m 2 , takes the form 



E ^ = ±TZ 1/2 , K = [E(r 2 + a 2 ) + L z a] 2 - A [m 2 r 2 + (L z + aE) 2 + Q] , (64) 

where A is the affine parameter of the trajectory. The conserved quantity Q takes the 
explicit form Q = L 2 — L 2 — a 2 (E 2 — m 2 ) cos 2 #00, where 9^ is the value of the azimuthal 
angle as the particle approaches the black hole from infinity, and L the total angular 
momentum of the particle. 

A particle approaching a rotating black hole from infinity, may do so by following 
a number of possible trajectories. Here, we will be interested in the case of a massless 
particle with its trajectory being either transverse (9^ = 7r/2) or parallel (9^ = 0, n) to 
the rotation axis. Starting with the first case, we notice that, for motion strictly on the 
equatorial plane, cos#oo = and L = L z . Then, Q = 0, and Eq. (1641) takes the form 

( E |)^[^ + »£ + ( ^-, 4 m 

In the above, we have defined L/E = b, where b > is the impact parameter of the 
particle. For the above equation to be consistent, its right-hand-side, or equivalently 
the expression inside the square brackets, must be positive-definite. Since the particle 
approaches the black hole from large r, we focus our attention on the radial regime outside 
the ergosphere, where the coefficient of b 2 , (a 2 — A), can be shown to be negative. Then, 
the constraint on the values of b takes the form 62 < b < b±, where 61,2 are the roots of 
the equation, following by setting the right-hand-side of Eq. (|65|) equal to zero. However, 
it may easily be seen that 6 2 < 0, therefore, the classically acceptable regime is defined 
by the constraint < b < b±. Particles with impact parameters in this regime can access 
all values of the radial coordinate, and thus reach the black-hole horizon, too, where they 
get absorbed. According to the geometrical optics argument, then, the closest distance 
the particle can get to the black hole without being captured is 



an + r n+1 J a 2 + r 2 — 



min(&i) = min . (66) 

1 r n+i _ ^ / v ' 

As a consistency check, we observe that, for a = 0, the above expression reduces to r c = 
min (r/ v/l — ^ht) > which leads directly to the result fl63l) for a non-rotating brane black- 
hole background derived in [HI [13]. By further setting n = 0, the purely 4-dimensional 
Schwarzschild case (401 EJUS] is also recovered, with r c = 3v^3t7i/2. 

For general n and a, an analytic expression for the minimum distance r c is difficult 
to find. Nevertheless, a simple numerical analysis may lead to the value of r c in units of 
rh, after using Eq. to eliminate the mass parameter \x from Eq. fl66|) . Then, through 
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a \ n 


1 


2 


3 


4 


5 


6 


u.u 


12.0 


9.6 


8.2 


1.6 


6.6 


6.2 


0.3 


17.9 


12.0 


in yl 

10.4 


9.0 


8.1 


7.4 


0.5 


23.5 


15.9 


12.6 


10.7 


9.5 


8.7 


0.7 


31.0 


20.0 


15.4 


12.9 


11.4 


10.3 


1.0 


46.0 


27.6 


20.6 


17.0 


14.8 


13.4 


1.5 


81.9 


44.3 


31.7 


25.6 


22.2 


19.9 


2.0 


131.6 


65.4 


45.3 


36.3 


31.3 


28.2 



a*\n 


1 


2 


3 


4 


5 


6 


0.0 


12.6 


9.6 


8.2 


7.3 


6.6 


6.2 


0.3 


10.0 


8.2 


7.3 


6.6 


6.2 


5.9 


0.5 


9.5 


8.1 


7.3 


6.8 


6.4 


6.1 


0.7 


9.5 


8.4 


7.7 


7.3 


6.9 


6.7 


1.0 


10.5 


9.6 


9.0 


8.6 


8.4 


8.1 


1.5 


13.4 


13.2 


12.7 


12.3 


12.1 


11.9 


2.0 


19.2 


18.5 


18.1 


17.8 


17.5 


17.4 



Table 2: Absorption cross-section a a bs (in units of r\) for particles moving in the equa- 
torial plane of the axially-symmetric brane black hole ([H]), for a > (upper sub-table) 
and a < (lower sub-table). 



the relation <7 a b s = nr 2 , the corresponding absorption cross-section may be found; its 
values, for a variety of n and a*, are displayed in Table 2. The two sub-tables correspond 
to the two possible orientations of the particle's angular momentum L: as it approaches 
the black hole from infinity moving in the equatorial plane, its angular momentum and 
the black-hole one can either be parallel {aL > 0) or anti-parallel (aL < 0). Here, we 
have assumed that L > always, and considered two different choices for the sign of the 
angular momentum parameter of the black hole, a > and a < 0, that correspond to 
the first and second sub-table of Table 2, respectively. For a < 0, the sign of the a/z-term 
in the numerator of Eq. fl6"ET) is reversed, a modification that leads to a lower value of r c 
and eventually of the cross-section. 

We now proceed to the case of a zero-mass particle coming from infinity in an orbit 
parallel to the black hole's rotation axis. This translates into cos 2 9^ — 1 and L z = 0. In 
that case, we find 

Z^-) =E 2 (r 2 + a 2 ) 2 -AL 2 . (67) 



Defining, as before, b = L/E, one may easily conclude that the above equation is again 
consistent only if 



The above leads to the minimum distance of the particle's approach to the black hole 
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a„\n 



1 



2 



3 



4 



5 



6 



0.0 
0.3 
0.5 
0.7 
1.0 
1.5 
2.0 



12.6 
13.7 
15.7 
18.7 
25.1 
40.9 
62.8 



9.6 

10.4 

11.6 

13.5 

17.2 

25.7 

36.5 



8.2 

8.7 

9.7 

11.2 

14.0 

20.5 

28.9 



7.3 

7.8 

8.6 

9.8 

12.3 

18.0 

25.3 



6.6 

7.1 

7.9 

9.0 

11.3 

16.5 

23.3 



6.2 

6.6 

7.3 

8.4 

10.5 

15.4 

22.0 



Table 3: Absorption cross-section <7 a b s , in units of r^, for particles moving parallel to the 
rotation axis of the axially-symmetric brane black hole ©. 

without being captured given by 



For a = 0, the above result also reduces to the Schwarzschild-like one ( 163|) . as expected, 
since, in the absence of rotation, all directions of the particle's orbit should give the same 
result. By using Eq. (1691 . the values of the corresponding absorption cross-section <r a b s , 
in units of r^, are given in Table 3. 

Let us now compare the various analytic values, obtained above by using the geo- 
metrical optics limit, with the numerical ones for the high-energy asymptotic value of 
the absorption cross-section. One should, of course, be careful when a direct comparison 
of these results is made: the values displayed in Tables 2 and 3 correspond to trajectories 
with a specific azimuthal angle, while the numerical values of Table 1 are actually inte- 
grated over all angles. Nevertheless, a comparison between the two sets of results could 
reveal, upon finding an agreement, which type of trajectories may be used to account for 
the value of the total cross-section more accurately. A direct comparison of the entries of 
Tables 1 and 2 shows that the total cross-section is smaller than the one corresponding 
to a trajectory lying on the equatorial plane with aL > 0, but larger than the one with 
aL < 0. On the other hand, by comparing the entries of Tables 1 and 3, we find that 
there is an almost perfect agreement between these two sets of results for low values of 
either a* or n. The same agreement can be pictorially seen in Fig. 0, where the partic- 
ular case n = 2 and a* = 0.3 is shown. We may thus conclude that particle trajectories 
running parallel to the rotation axis of the black hole lead to an absorption cross-section 
whose value is virtually identical to the total one. As either n or a* increases further, 
the two sets deviate; in this parameter regime, the contribution of all possible particle 
trajectories needs to be more carefully taken into account before the value of the total 
absorption cross-section can be justified. 




(69) 
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5 Conclusions 



The emission of Hawking radiation, i.e. the evaporation of a black hole via the emission 
of elementary particles, takes place during the spin-down and Schwarzs child phase of 
its life. Although the emission during the Schwarzschild phase of a higher- dimensional 
black hole was studied, both analytically and numerically, quite early, the complexity 
of the gravitational background around a similar, but rotating, black hole delayed the 
study of the spin-down phase. During the last few years, numerical studies have derived 
results for the various spectra characterising the emission of elementary particles on 
the brane by a higher-dimensional rotating black hole - the most phenomenologically 
interesting emission channel for the brane-localised observers. Nevertheless, no analytical 
studies have been performed and no analytic expressions for the fundamental quantities 
governing the emission of Hawking radiation, such as the absorption probability, have 
ever been derived for an arbitrary value of the number of extra dimensions n. In this 
work, we have duly performed this task, and studied in detail the properties of the 
absorption probability and absorption cross-section for scalar fields emitted on the brane 
by the (4 + n)-dimensional axially-symmetric black hole. 

As the complexity of the equation of motion, describing the propagation of a scalar 
field in the axially-symmetric brane background, forbids the derivation of a general so- 
lution for a particle with arbitrary frequency, we were forced to focus our analysis on 
two particular energy regimes: the low-energy one and the high-energy one. The low- 
energy regime was studied in Section 3, where an analytic solution for the radial part 
of the scalar-field wave-function was derived. This involved matching the near-horizon 
and far-field asymptotic solutions in an intermediate regime, and allowed calculation of 
the absorption probability. Our analytic results, formally valid only for low values of the 
energy parameter urh and angular momentum parameter a*, were compared with the 
exact numerical results, that were also reproduced during our analysis in an attempt to 
check the range of validity of our approximations. The two sets of results were found to be 
in excellent agreement in the low-energy regime, as expected. In addition, even for mod- 
erately large values of ojr h and a*, the agreement on both qualitative and quantitative 
levels still persisted. 

The properties of the absorption probability in the low-energy regime, as these follow 
from our analytic results, were then studied in detail. Its dependence on the angular mo- 
mentum numbers (/, m) was investigated first, with our analysis revealing the dominance 
of the lowest partial mode I = 0, and the one with m — — \l\, for fixed I. The absorption 
probability was also found to strongly depend on the spacetime topological parameters, 
namely the angular momentum parameter a* of the black hole and the number of extra 
dimensions n: the non-superradiant modes with m < 0, were shown to be enhanced 
with both a* and n, while the superradiant ones, with m > 0, were, on the contrary, 
suppressed - however, outside the superradiant regime, the latter modes were enhanced 
with n. Our analytic expression for the absorption probability, valid in the low-energy 
regime, was then expanded in the limit u — > 0, and the constant asymptotic value of the 
absorption cross-section was derived. The analytic value was shown to accurately repro- 
duce, for small a*, the exact numerical one, that was equal to the area of the horizon of 
the projected-on-the-brane axially-symmetric black hole, 47r(a 2 + r%). 
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We subsequently turned our attention to the study of the scalar equation of motion 
in the high-energy regime. By using an analogous approximate method, and for the first 
time in the literature, an analytic solution was derived that perfectly matched the exact 
numerical one for large enough values of the energy parameter lor^. The value of urh, 
beyond which the two solutions completely coincide, was shown to be strongly mode- 
dependent: modes with small I allowed the two solutions to match fairly quickly, while, 
modes with larger / had the matching taking place at an increasing value of the energy 
parameter. By employing the exact numerical solution for the absorption probability, 
valid at all energy regimes, we were able to determine the constant asymptotic value of 
the total absorption cross-section at the high-energy regime, and its dependence on n and 
a*. Similarly to the spherically-symmetric case, this asymptotic value is suppressed as 
the number of extra dimensions increases. On the other hand, an increase in the angular 
momentum of the black hole causes an enhancement in the high-energy asymptotic value 
of cr a bs- A detailed analysis, based on the geometrical optics limit, revealed that the 
asymptotic value of the absorption cross-section in the high-energy limit is accurately 
reproduced by considering particle trajectories approaching the black hole from infinity 
and running parallel to the rotation axis of the black hole. 

The analytic results, supplemented by exact numerical ones, derived in this work, on 
the behaviour of the absorption probability and cross-section for scalar particles propa- 
gating in an axially-symmetric brane black-hole background, smoothly complement the 
sole previous analytic study of the 5-dimensional case [32] , as well as the numerical stud- 
ies of the Hawking radiation spectrum of [231 USE [27]. Given the excellent agreement 
between our analytic solutions and the exact numerical ones, in the low and high-energy 
regime, their use to derive the corresponding emission rates would have led to results 
identical to the ones already presented in the works cited above. For that reason, we 
have refrained from performing this task here, and refer the interested reader to those 
works. Instead, in this manuscript, we have focused our attention on the derivation of 
closed-form expressions and study of the properties of the absorption probability and 
cross-section, that carry a significant amount of information on particle properties as 
well as on properties of the spacetime. Apart from their obvious theoretical interest, the 
above results may have a phenomenological one too: both our analytic solutions could be 
reliably used, in place of the exact numerical ones, to derive the energy emission rates in 
the low and high-energy regimes; these could then be used for the interpretation of any 
observable effects coming from an evaporating black hole produced in a ground-based 
accelerator and centered in these two frequency regimes. In the high-energy regime, in 
particular, where constraints on the available running time may put limits on the deriva- 
tion of numerical data, our analytic solution may prove useful in removing the need for 
numerical integration. 
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